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Summary 

A computational study has been conducted in which the CFL3D Navier-Stokes solver coupled with an 
algebraic and an one-equation nonequilibrium turbulence model has been used to predict the flow over a 
65* delta wing at transonic conditions for Reynolds numbers ranging from 6x1 0 6 to 120x10 6 based on 
mean aerodynamic chord. Solutions obtained indicated that the computational method when used with 
the one-equation turbulence model predicts results that compare well with experiment for attached flow 
conditions. Comparisons with experimental pressures at separated conditions show that the 
computational method, even though primary flowfield features are predicted well, does not predict 
secondary flow features. 

Introduction 

The flow over low-aspect ratio wings exhibit some very interesting flow features. The flow around the 
forebody, strake and highly swept leading-edges of fighter aircraft such as the F-18 1 and the F-106B 2 sep- 
arates at moderate to high angles-of-attack. The resulting free shear layer rolls up in a spiral fashion to 
form the leading-edge vortex. Locally high velocities are induced by the vortex in the flow field that in turn 
creates a low pressure region near the vehicle surface. This suction peak located on the wing surface 
under the vortex core results in additional lift that is known as “vortex lift”. Vortices from the forebody, 
strake and the wing leading-edge create a very complex flow field for fighter type configurations while for a 
sharp edge delta wing planform only a single primary vortex forms. However, this flow field is not as sim- 
ple as it may seem. An adverse lateral pressure gradient exists outboard of the minimum pressure region 
created by the primary vortex. This gradient results in a secondary separation creating another vortex, 
weaker than the primary, turning in the opposite sense of the first. When the secondary vortex is relatively 
strong (high angles-of-attack), an adverse lateral pressure gradient now exists inboard of the secondary 



vortex resulting in a third separation and a tertiary vortex. (Hoeijmakers gives an excellent overview of the 
related physics of vortical flowfields in references 3 and 4.) 

Aerodynamic benefits can be obtained from vortical flow when the wing vortex is somewhat 
controlled 5 ' 6,7 ’ 8 . In the applications presented in these references, a vortex is induced by fixing the separa- 
tion point along the sharp leading-edge of a leading-edge flap or “vortex flap”. The vortex forms on this 
flap resulting in a beneficial force vector that can be resolved into thrust and lift components. A naturally 
occurring, uncontrolled, vortex may result in adverse effects due to asymmetric separation along the fore- 
body, strake or wing leading-edge 9 . These effects can cause buffet or a loss of stability and control, and 
depend on how early or late the flow separates. Flow separation around a rounded leading-edge is a flow 
phenomenon that is very dependent on Reynolds number and, since boundary layer transitional data is not 
a linear function of Reynolds number, this information can not be accurately predicted by standard (low 
Reynolds number) wind tunnel testing. The capability to accurately predict this separation is gaining signif- 
icance in the design of some current configurations such as the High Speed Civil Transport (HSCT) which 
will consist of a highly swept wing with a round leading-edge. 

Numerous experimental as well as numerical studies have been conducted to gain an understanding 
of this flow field. Quite a few of these studies have been done using a delta wing planform since, as stated 
above, this simple geometric shape can produce the desired complex vortical flow field. (See references 
10 and 11 for experimental work and references 12-26 for computational studies.) The experimental 
studies were all conducted at sub-scale Reynolds number conditions so the numerical studies concen- 
trated on geometries with sharp leading-edges. These studies, most either Euler or laminar Navier-Stokes 
solutions, predicted the primary vortex trajectory and strength well while not capturing the secondary vor- 
tex. This is expected since correctly capturing the secondary vortex depends upon capturing the second- 
ary separation, a viscous phenomenon, while for a sharp leading-edge the primary vortex forms as a result 
of the geometrical discontinuity at the wing leading-edge. A few researchers have conducted numerical 
studies on rounded leading-edge geometries. Hsu and Liu 24 and Hsu, Hartwich and Liu 25 showed compu- 
tations that compared well with experiment for laminar, incompressible conditions. Rizzi and Purcell 13 and 
Muller and Rizzi 14 showed laminar Navier-Stokes solutions that compared well with experiment for predic- 
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tion of the primary vortex at transonic conditions but failed to completely capture the secondary vortex. 
Hilgenstock 12 and Hilgenstock and Vollmers 21 made fully turbulent Navier-Stokes calculations using an 
algebraic turbulence model at transonic and moderately low Reynolds number. Their prediction of the pri- 
mary vortex compared well with experiment while the secondary separation was not predicted. 

Recently an experimental investigation of the flow over a 65 delta wing was conducted in the National 
Transonic Facility (NTF) at the NASA Langley Research Center 27 . Pressure data is available from this 
wind tunnel test at transonic conditions for Reynolds number of up 120 million, based on mean aerody- 
namic chord. This data was obtained for a variety of leading-edge radii ranging from sharp to fairly blunt. 
This test was conducted to establish a high Reynolds number database so that the effect of leading-edge 
shape on leading-edge flow separation could be analyzed both numerically and experimentally. 

This paper will document a study in which the CFL3D Navier-Stokes code 28 was evaluated for the 
prediction of the flow over the NTF delta wing with a medium and a round radius leading-edge. These 
calculations were made for low to high Reynolds number at transonic conditions 

span 
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freestream Mach number 
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grid coordinate in the axial direction 
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z 


normal dimension, inches 


y + law-of-the-wall coordinate 

a angle of attack, degrees 

A sweep angle, degrees 


Subscripts: 

LE Leading-edge 

MAC Mean Aerodynamic Chord 

r root 

w wing 

oo denotes conditions at infinity 

Experimental Study 

An experimental investigation was conducted in the National Transonic Facility (NTF) at the NASA 
Langley Research Center in which pressure coefficients were obtained for a delta wing model with a series 
of leading-edge radii. The NTF delta wing model is an analytically defined flat plate delta wing with a lead- 
ing-edge sweep of 65° , 12 inch semi-span and 25.734 inch root chord (figure 1). This model was designed 
such that the leading-edge portion could be interchanged. Data was obtained for four leading-edge 
shapes in this test. These shapes varied from sharp to blunt and are analytically defined. Figure 2 shows 
the relative radii of the four leading-edge shapes. Data were obtained for each leading-edge radius for 
Reynolds number of 6x1 0 6 to 120x1 0 6 , based on mean aerodynamic chord at freestream Mach numbers 
of 0.40 to 0.90. All experimental data was obtained at free-transition conditions. Angle-of-attack ranged 
from O' to 26*. The wind tunnel test data is reported in Reference 27. 

For the present study, analyses were preformed for the wing with the medium and round leading-edge 
radius. The leading-edge shape is defined by 
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z 


a.7x + bx + cx 2 + dx 3 


(1) 


where a = 0.34897572 , b = 0.29224873 , c = -0.28382228 , and d = 0.06227092 for the medium 
radius leading-edge and a = 0.4935262 , b = 0.08004809 , c = -0.19719469 .and d = 0.04632388 

for the round leading-edge radius. The radius for each of these leading-edges non-dimensionalized by the 
mean aerodynamic chord is 0.15% and 0.30% for the medium and round leading-edges, respectively. The 
trailing-edge closure can also be defined analytically by 

z = bx + cx 2 + dx 3 (2) 

where b = 0.510024 , c = -0.19819009 , and d = 0.02567154 and is the same for each of the 
leading-edge shapes. The sting cover faring is defined by 

z = aVx + bx + cx 2 + dx 3 (3) 

where a = 0.30946082 , b = 0.33279823 , c = -0.04163681 .and d = 0.00150729 

Computational Method 

The computational method used in this study is the CFL3D, version three, Navier-Stokes code 28 . 

This code, developed at NASA Langley, solves the thin-layer approximations to the three-dimensional, 
time-dependent, compressible Navier-Stokes equations. A second-order accurate upwind-biased spatial 
differencing scheme is utilized for solving these equations. Convergence acceleration is obtained using 
the full approximation scheme multigrid algorithm and Roe’s flux-difference splitting scheme to determine 
fluxes on cell faces. The Baldwin-Lomax algebraic turbulence model 28 and the Spalart-Allmaras 
one-equation turbulence model 30 were used to solve for the turbulent viscosity and are incorporated into 
the laminar Navier-Stokes equations. The details of the numerical method as well as of each turbulence 
model are presented in the respective references. 

Grid Generation 

The grid generated about the delta wing geometry was created using the GRIDGEN grid generation 
package 31 . A nine block grid was created about the experimental geometry including the sting and the 
sting/wing fairing (see figure 3 for schematic of the block layout). The resulting mesh had nearly 900,000 
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grid points. Clustering normal to all surfaces was set to 1x10 6 c r This spacing was selected to ensure y + 
values at first grid line off the surface of 1 , based upon flat plate boundary layer theory, for the 
Rn = 120x1 0 6 case. All outer boundaries were set 20c r away from the geometry surfaces. Details of the 
surface grid as well as a crossflow grid plane are shown in figures 4 and 5. This gird was developed 
utilizing the arbitrary patching logic in CFL3D as is evident in these figures. 

Results 

Grid Resolution 

The effect of mesh density on the computed pressure coefficients was investigated by analyzing 
solutions on suscessively coarser grids. Each coarser grid was developed by eliminating every other point 
in each of the three coordinate directions. The resulting pressure distributions revealed that, while the 
upper surface results varied, the calculated lower surface pressure coefficients were insensitive to mesh 
density (figure 6). Proceeding on this information, a new volume grid was developed in which grid points in 
the streamwise direction were removed from several of the lower surface blocks and were added to the 
leading-edge block in the radial direction. The number of grid points around the leading-edge was 
increased from 49 to 73 in an attempt to better resolve the separation line. The total mesh size was kept 
fairly constant by deleting the grid points from the lower surface. The resulting pressure distribution 
exhibited in figures 7. a - 7.c shows that refinement of the leading-edge in the radial direction has little effect 
on the calculations. Since solutions using this grid did not utilize significantly more computer resources 
while yielding slightly better results it was selected to replace the original baseline grid (see table 1 for a 
listing of the original and refined grid dimensions). A successive coarsening of this grid yielded solutions 
that showed the upper surface pressure were not grid converged (figure 8). Refinement of the grid in the 
streamwise direction in the leading-edge block by doubling the number of grid points in the i direction and 
refinement of block 1 ,of block 7 , and of block 5 by doubling the number of j grid points in each of these 
blocks showed little to no effect on the calculated upper surface pressures (figure 9). This tends to indicate 
that the new baseline grid is grid converged in the radial and streamwise directions (see figure 10 for the 
convergence histories plotted for each of these grid refinement solutions). Solution analysis shows that 
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the maximum y + value utilizing the new baseline grid at the conditions that this grid resolution study has 
been made (M = 0.85, a = 16.37", and Rn = 6x1 0 6 ) were not greater than 1.60 with an average of 0.30 
and standard deviation of 0.18 indicating that the grid is sufficiently clustered to the surface in the normal 
direction. (A listing of the maximum and average with the corresponding standard deviation of y + for all 
solutions obtained utilizing this refined grid is presented in table 2.) 

Turbulence Model Comparison 

Comparisons of pressure coefficient predictions obtained using the Baldwin-Lomax turbulence model 
with those obtained using the Spalart-Allmaras turbulence model are presented in figures 11 .a - 11 c. 
Pressures predicted by either of the models compared very well with experiment for the attached flow case 
(a = 7.15'). There is slight if any difference between the predicted pressure coefficients except near the 
trailing-edge (figure 11a). Here the solution obtained using the Spalart-Allmaras turbulence model 
exhibits pressure distributions that more closely follow the experimental pressure distribution. Comparison 
of predicted pressures using both turbulence models with experiment is nearly as good for the case with 
slight flow separation at the leading-edge of the wing (a = 12.3*). Pressures obtained using both models 
show evidence of a delayed leading-edge separation (figure ll.b, x/c r = 0.20). The location of the 
experimental primary vortex pressure peak is predicted equally as well using either turbulence model 
however the primary vortex suction peak is predicted better using the Baldwin-Lomax turbulence model 
(sixty and eighty per cent root chord). The onset of leading-edge separation is also not predicted for the 
fully separated case (a = 16.37*) by either of the turbulence models as shown by the prediction of the 
formation of the pressure suction peak aft of experiment as indicated by the delay in the formation of the 
computed pressure peak (figure ll.c, x/c r = 0.20). The solutions obtained using the Spalart-Allmaras 
turbulence model better predicts the magnitude of the suction peak although the crossflow shock wave 
(defined by (Ac p ) / (A (y/ (b/2) ) ) > 8 outboard of the primary suction peak) evident in the experimental 
pressures aft of 60% root chord (approximately 65% semi-span) is not predicted. Pressure coefficients 
predicted with the Baldwin-Lomax turbulence model compared very poorly with experiment at 80% root 
chord. Investigation of this solution further showed that the pressure coefficients predicted aft of the 80% 
root chord location were oscillating (figure 1 2). Since this behavior was only observed with use of the 
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Baldwin-Lomax turbulence model the Spalart-Allmaras turbulence model was selected for the remainder of 
this study. 

Solution Analysis 

Predicted pressure coefficients compared with experiment are shown in figures 11 .a - 11 ,c and in 
figures 13. a - 15.c for the medium leading-edge geometry and in figures 16.a - 19.c for the round 
leading-edge geometry. The pressures predicted using the Spalart-Allmaras turbulence model compare 
well with experiment for the attached flow case (low angles-of-attack) for all Reynolds numbers computed. 
As angle-of-attack is increased the flow separates near the apex of the wing in the experiment resulting in 
a pronounced suction peak indicating a leading-edge vortex (e.g., figure 11. b, x7c r = 0.40, y/(b/2) = 0.80). 
The vortex is predicted in the calculations but the location of the formation of the suction peak is predicted 
aft of the experimental location. A secondary separation is evident in the experimental pressures as 
indicated by the formation of an additional suction peak outboard of the main suction peak at the highest 
angle-of-attack analyzed at eighty per cent root chord, eighty per cent semi-span for the medium radius 
leading-edge geometry (see figures 11c, 13. c, 14. c and 15.c) and at x/c r = 0.60, y/(b/2) = 0.80 for the 
round radius leading-edge geometry (see figures 16.c, 17.c, 18.c and 19.c). This separation is not 
predicted using the computational method as shown by the surface streamlines plotted in figures 20.a - 
21. d. In these plots, a separation line is indicated by the streamlines coalescing to a line while an 
attachment line is a line from which the streamlines are flowing out. As can be seen in the streamline plots 
the computational method is predicting only a single separation and corresponding attachment line. 

Experimental pressure distributions show little effect of Reynolds number for both leading-edge 
shapes, particularly at low angles-of-attack where the flow is attached. The Navier-Stokes code generally 
predicts the experimental trend of Reynolds number on pressure distributions (figures 22.a - 23.f) with the 
exception of the Rn = 48x10 6 for the medium leading-edge geometry with a = 12‘. Here, experiment is 
showing an attachment line (the increase in pressure inboard of the suction peak) that is much nearer to 
the leading-edge than for the other Reynolds numbers analyzed (figure 22.c, x/c r = 0.40) as well as a lower 
suction peak for the Rn = 48x1 0 6 data (figure 22. c, x/c r = 0.60). The computational method predicts 
neither flow feature for the Rn = 48x1 0 6 case (figure 22. d). Further analysis of the experimental data 
presented in Reference 27 for the medium radius leading-edge geometry reveals that the trend exhibited 
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in the experimental pressures for a = 12.03” is entirely a Reynolds number effect, an effect that the 
Navier-Stokes code utilized in this study did not predict (compare figures 22.c and 22. d). Also evident in 
the experimental pressure distributions for Rn = 48x1 0 6 is the formation of a dual leading-edge vortex as 
indicated by the dual pressure peaks of near equal magnitude (figure 21 c, x/c r = 0.40, y/(b/2) > 0.80). 
These pressure peaks are not as prominent in the pressure distributions for Rn = 6x10 or 
Rn = 84x1 0 6 . The dual vortices are not predicted by the computational method. A crossflow shock 
wave (defined by (Ac p ) / (A (y/ (b/2) ) ) > 8 outboard of the primary pressure peak) is evident in the 
experimental pressures (e.g., figure 16.c, x/c r = 0.60, y/(b/2) = 0.65). This flow phenomena is not 
predicted by the Navier-Stokes code. Lower surface pressure coefficients were predicted very well at all 
angles-of-attack and Reynolds numbers indicating that the angle-of-attack was correctly modeled in the 
calculations. 

Comparison with data derived from an empirical relationship presented in reference 32 shows that the 
predicted sweep and angle-of-attack of the vortex core agree very well with experimentally derived values 
(table 3). Predicted values were obtained by locating the vortex in the flowfield by identifying the minimum 
total pressure above the boundary layer over the wing at stations along the chord of the wing. 

Examination of the minimum pressure coefficient at x/c r = 0.60 in the footpath of the vortex core (figures 24 
and 25) shows that the experimental level is generally under-predicted. The experimental trends are not 
predicted for either of the leading-edge radii analyzed with the predicted data showing very little effect of 
Reynolds number. (Note that the Reynolds number effect on the medium radius leading-edge for values 
between Rn = 6x1 0 6 and Rn = 84x1 0 6 , as discussed above, is again evident in the experimental data 
(figure 24, a = 12”)). 

Concluding Remarks 

This study has “indicated" several conclusions regarding the use of Navier-Stokes methods with 
turbulence modeling to predict the flow over a highly swept wing. Among these is that, at attached flow 
conditions, the predictions of the Navier-Stokes method employed agree very well with experiment. 
However, when the flow is fully separated comparisons with experimental data are not as good. This is a 
result of several factors, individually or combined. The first of these is that the solutions were computed for 
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fully turbulent flow conditions while experimental data was obtained for free transition conditions. Further 
complicating matters, this data was obtained in the National Transonic Facility in which there are no means 
presently available to obtain flow visualization data that could indicate where transition was occurring in 
the experiment. This lack of experimental information could be offset with the development of some 
affordable means to calculate transition locations that can be coupled with current Navier-Stokes solvers. 

Comparison of experimental data with Navier-Stokes solutions are also degraded due to the 
computational method not predicting any of the secondary flow features apparent in the experimental 
pressure distributions. The location of the primary separation and attachment are predicted well although 
the streamwise location of leading-edge separation is not. The location of the primary vortex suction peak 
is predicted well thus indicating that the computational method is modeling the overall flowfield physics 
correctly, which is as expected. The Navier-Stokes code does not predict the secondary flow features 
such as the dual leading-edge vortices and the crossflow shock wave. These features could be captured 
in several ways. One of these is through grid enhancement in the region of the vortex. Increasing grid 
density in this region should allow for better resolution of the vortical flowfield. However, in a design 
environment these flow features are often not known a priori. Therefore for a computational method to be 
a useful design tool some means of adapting the grid to the vortex is required. Another method by which 
the vortical flow could be better resolved is through better turbulence modeling techniques. The current 
nonequilibrium turbulence model utilized in this study predicts primary separation and attachment well as 
well as lower surface pressures for all cases analyzed. This model has difficulties in predicting secondary 
flow features for the geometry analyzed. 

This study has shown that the CFL3D Navier-Stokes code coupled with the Spalart-Allmaras 
one-equation nonequilibrium turbulence model can predict the primary flow features of a highly swept wing 
in a transonic flowfield as well as first order Reynolds number effects. However, absolute pressure 
distribution levels are not predicted, a consequence of the theory not predicting any of the secondary flow 
features. With this in mind, this computational method can be used to predict first order effects on the 
flowfield. 
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Figure 3. Schematic of block layout. 
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Figure 4. Head-on view of computational grid. 



Figure 5. Close-up of crossflow computational grid plane, medium leading-edge radius. 
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Figure 6. Lower surface pressure coefficients show no dependence on mesh density: Medium 
leading-edge radius, Spalart-Allmaras turbulence model, = 0.85, a = 16.37", Rn = 6x1 0 6 . 
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Figure 7. a. Refinement of leading-edge block shows little effect on pressure coefficients: Medium 
leading-edge radius, Spalart-Allmaras turbulence model, M„ = 0.85, a = 7.15 , Rn = 6x10 . 
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Figure 7.b. Refinement of leading-edge block shows little effect on pressure coefficients: Medium 
leading-edge radius, Spalart-Allmaras turbulence model, M m = 0.85, a = 12.3*, Rn = 6x1 0 6 . 
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Figure 7.c. Refinement of leading-edge block shows little effect on pressure coefficients: Medium 
leading-edge radius, Spalart-Allmaras turbulence model, = 0.85, a = 16.37‘, Rn = 6x10 . 
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Figure 8. Refined grid upper surface pressures are not grid converged: Medium leading-edge radius, 
Spalart-Allmaras turbulence model, M«, = 0.85, a = 16.37’, Rn = 6x1 0 6 . 
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Figure 9. Upper surface pressure distributions show little effect due to upper surface grid refinement: 
Medium leading-edge radius, Spalart-Allmaras turbulence model, M„ = 0.85, a = 16.37", Rn = 6x10 . 
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Figure 11. a. Turbulence model comparisons: 

Rn 
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Medium leading-edge radius, M m = 0.85, a = 7.15*, 
= 6x1 0 6 . 
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Figure ll.b. Turbulence model comparisons: Medium leading-edge radius, = 0.85, a = 12.3*, 

Rn = 6x1 0 6 . 
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Figure 11 .c. Turbulence model comparisons: Medium leading-edge radius, M„ = 0.85, a = 16.37°, 

Rn = 6x1 0 6 . 
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Figure 12. Baldwin-Lomax Turbulence model predicts oscillating pressures: Medium leading-edge 

radius, M„ = 0.85, a = 16.37°, Rn = 6x1 0 6 . 
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Figure 13.b. Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M OT = 0.85, a = 12.03', Rn = 48x10 . 
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Figure 14.a. Medium leading-edge radius pressure coefficients, Soalart-Allmaras turbulence model, 

M«, = 0.85, (X = 6.99\ Rn = 84x10 . 
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Fiaure 14 b Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

= 0.85, 0 C = 12.31', Rn = 84x10 . 
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Figure 14.c. Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

= 0.85, a = 16.59*, Rn = 84x10 . 
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Fiqure 15. a. Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M«, = 0.85, a = 7.03*. Rn = 120x10° . 
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Figure 15.b. Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

= 0.85, a = 12.44‘, Rn = 120x10 . 
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Figure 15.c. Medium leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M„ = 0.85, a = 14.63*, Rn = 120x10 . 
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Figure 16.a. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M„ = 0.85, CL = 7.14‘, Rn = 6xKJ* 
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Figure 16 .c. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M„ = 0.85, a = 16.37”, Rn = 6x10° . 
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Fiqure 17 a. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M =0.85, a = 6.98', Rn = 48x10°. 
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Figure 17. c. Round leading-edge radius pressure coefficients, = 0.85, a = 16.44*, Rn = 48x10 6 . 
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Figure 18 b Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

ML = 0.85, a = 12.09", Rn = 84x10 . 
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Figure 1 8.c. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

= 0.85, a = 1 6.38*, Rn = 84x1 0 . 
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Figure 19.a. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model 

= 0.85, a = 7.03\ Rn = 120x10 . 
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Figure 19.b. Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

M m = 0.85, a = 12.29', Rn = 120x10 . 
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Fiqure 19c Round leading-edge radius pressure coefficients, Spalart-Allmaras turbulence model, 

= 0.85, a = 16.59", Rn = 120x10 . 


47 













Figure 20.a. Computational surface streamlines, Spalart-Allmaras turbulence model, medium 
leading-edge radius, M„ = 0.85, a = 16.37', Rn = 6x10 . 



Figure 20.b. Computational surface streamlines, Spalart-Allmaras turbulence model, medium 
leading-edge radius, M m = 0.85, a = 16.25', Rn = 48x10 . 
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Figure 20. c. Computational surface streamlines, Spalart-Allmaras turbulence model, medium 
leading-edge radius, = 0.85, a = 16.59*, Rn = 84x10 . 



Figure 20. d. Computational surface streamlines, Spalart-Allmaras turbulence model, medium 
leading-edge radius, = 0.85, a = 14.63’, Rn = 120x10° . 
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Figure 21 .a. Computational surface streamlines, Spalart-Allmaras turbulaice model, round 
leading-edge radius, = 0.85, a = 16.37*, Rn = 6x10 . 



Figure 21 .b. Computational surface streamlines, Spalart-Allmaras turbulence model, round 
leading-edge radius, = 0.85, a = 16.44*, Rn = 48x10 ’ . 
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Figure 21. c. Computational surface streamlines, Spalart-Allmaras turbulence model, round 
leading-edge radius, = 0.85, a = 16.38', Rn = 84x10 . 



Figure 21. d. Computational surface streamlines, Spalart-Allmaras turbulence model, round 
leading-edge radius, = 0.85, a = 16.59*, Rn = 120x10 ’ . 
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Figure 22.a. Reynolds number effect on pressure coefficients, medium leading-edge radius, 

= 0.85, a = 7*. experiment. 
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Figure 22.b. Predicted Reynolds number effect on pressure coefficients, medium leading-edge radius, 
M m = 0.85, a = 7\ CFL3D, Spalart-Allmaras turbulence model. 
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Figure 22.c. Reynolds number effect on pressure coefficients, medium leading-edge radius, 

= 0.85, a = 12', experiment. 
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Figure 22. e. Reynolds number effect on pressure coefficients, medium leading-edge radius, 

M„ = 0.85, a = 16', experiment. 
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Figure 22.f. Predicted Reynolds number effect on pressure coefficients, medium leading-edge radius, 
= 0.85, a = 16°, CFL3D, Spalart-Allmaras turbulence model. 
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Figure 23. c. Reynolds number effect on pressure coefficients, round leading-edge radius, M 

a = 12", experiment. 

















Figure 23.d. Predicted Reynolds number effect on pressure coefficients, round leading-edge radius, 
= 0.85, a = 12”, CFL3D, Spalart-Allmaras turbulence model. 
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Figure 23.f. Predicted Reynolds number effect on pressure coefficients, round leading-edge radius, 
= 0.85, a = 16”, CFL3D, Spa art-Allmaras turbulence model. 
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Figure 24. Reynolds Number effect on minimum suction pressure for medium leading-edge radius 

geometry, = 0.85, x/c r = 0.60. 



Figure 25. Reynolds Number effect on minimum suction pressure for the round leading-edge radius 

geometry, = 0.85, x/c r = 0.60. 
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Block 

Original Grid 

Refined Grid 

! max 

Jmax 

^max 

J max 

Jmax 

^max 

1 

49 

13 

65 

49 

13 

65 

2 

49 

49 

65 

49 

73 

65 

3 

49 

13 

65 

33 

13 

65 

4 

49 

65 

65 

49 

65 

65 

5 

25 

49 

65 

25 

49 

65 

6 

49 

33 

65 

49 

33 

65 

7 

49 

25 

65 

49 

25 

65 

8 

49 

25 

65 

33 

25 

65 

9 

49 

33 

65 

33 

33 

65 


Table 1 : Computational grid dimensions. 
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Leading-Edge 

Radius 

Reynolds 

Number 

(xIO' 6 ) 

«w 

a core 
(ref. 32) 

a core 

(CFL3D) 

'^core 

(ref. 32) 

■\:ore 

(CFL3D) 

Medium 

6 

12.3’ 

4.07* 

3.96* 

73.8* 

73.81’ 

Medium 

6 

16.37’ 

5.09* 

5.10’ 

73.8* 

72.61’ 

Medium 

48 

12.03* 

4.01* 

4.23* 

73.8* 

74.02* 

Medium 

48 

16.25’ 

5.06’ 

5.10’ 

73.8* 

72.60’ 

Medium 

84 

12.03’ 

4.01* 

3.96’ 

73.8’ 

73.00’ 

Medium 

84 

16.25’ 

5.06’ 

5.10’ 

73.8* 

72.60* 

Medium 

120 

12.44’ 

4.11* 

3.77’ 

73.8’ 

73.12’ 

Medium 

120 

14.63* 

4.66* 

3.90* 

73.8* 

72.87* 

Round 

6 

12.28* 

4.07* 

4.05* 

73.8* 

73.41* 

Round 

6 

16.37* 

5.09* 

5.00* 

73.8* 

73.11’ 

Round 

48 

12.24* 

4.06* 

4.42* 

73.8’ 

74.43* 

Round 

48 

16.44* 

5.11* 

5.00’ 

73.8’ 

73.11’ 

Round 

84 

12.09’ 

4.02’ 

3.30* 

73.8’ 

73.80* 

Round 

84 

16.38’ 

5.09’ 

5.00’ 

73.8’ 

73.11’ 

Round 

120 

12.29’ 

4.07* 

3.17* 

73.8’ 

73.90* 

Round 

120 

16.59’ 

5.15* 

5.00* 

73.8* 

73.1 1* 

1 


a core = a w " ( (01 + 0 01 x A LE ) oi w - 1.0) , deg. (Equation 3 from Reference 32) 
A CO re = 33.5 + 0.62A le , deg. (Equation 2 from Reference 32) 



Table 3: Predicted vortex core trace compares well with experiment. 
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